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ABSTRACT 


We consider electron acceleration by obliquely propagating fast mode waves 
in magnetically dominated accretion disk coronae. For low coronal plasma 
densities, acceleration can exceed Coulomb drag at lower energies and energize 
electrons out of the thermal background, resulting in a nonthermal tail. The 
extent of this tail is determined by the balance between acceleration and 
radiative cooling via inverse Compton scattering and synchrotron emission, and 
usually goes out to tens of MeV. This will have direct applications for explaining 
the gamma-rays from several galactic black hole candidates, such as Cyg X-1 
and GRO J0422, which show 0.5-5 MeV emissions in excess over what most 
thermal models predict. Detailed time evolutions of the particle distributions 
and wave spectra are also presented. 

Subject headings: acceleration of particles — accretion, accretion disks — 
gamma rays: theory — waves 
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1. Introduction 


Thermal Comptonization models have had much success in explaining the hard X-ray 
spectra from galactic black hole candidates (GBHCs) with the plasma temperature ~ 50-100 
keV and Thomson depths of a few (e.g., [Shapiro, Lightman, &: Eardly 1976| ; punyaev "S 
Titarchuk 1980| ; [Harmon et ah 1994] ; [Liang 199^ . However, the most sensitive observations 
of GBHCs to date in the 0.5-5 MeV range by Compton Gamma-Ray Observatory have 
clearly revealed that persistent gamma rays (> 1 MeV) are being produced in some GBHCs, 
notably Cyg X-1 ([McConnell et ah 199ti[ ; [Phlips et al. 199^ ; [Ling et al. 1996[) and GRO 
J0422 ( yan Dijk et ah 1995[) . These gamma-ray emissions are very difficult to accommodate 
by the pure thermal models (e.g., [Sunyaev fc Titarchuk 1980| ; [Titarchuk 1994[) , strongly 
suggesting the need for modihcation ([Skibo fc Dermer 1995[) or to incorporate some 
nonthermal processes. 

Recently, Li, Kusunose, & Liang (1996) have proposed a gyroresonant stochastic 
electron acceleration model to account for the MeV emissions from GBHCs (see also 
Dermer, Miller, & Li 1996). In that model, they showed that high frequency whistlers 
can accelerate electrons directly from the thermal background, after which Alfven waves 
would continue the acceleration to higher energies. The result is an electron distribution 
with a hard non-Maxwellian tail. Compton scatterings from both thermal and nonthermal 
electrons produce a broad band X- to gamma-ray spectrum, in agreement with that 
observed from both Cyg X-1 and GRO J0422. One uncertainty associated with that 
model is the generation of whistlers, which can arise from either a cascade of wave energy 
from lower frequencies or a microinstability (see e.g., [Gary 19M[ ). While the nature of 
wave cascading is fairly well known in the MHD regime (e.g., [Verma 1994[ [Roberts'^ 


Wiltberger 1995 ), it has not been investigated at higher frequencies. Producing waves that 
can gyroresonate with an electron of a given energy through a resonant microinstability 
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requires (among other things) an anisotropic electron distribution containing electrons of 
that energy. Nonresonant instabilities may also produce the needed waves, but still require 
specihc anisotropic distributions. 


In this Letter, we consider electron acceleration by MHD fast mode waves that likely 
exist in accretion disk coronae. These low-frequency waves can accelerate electrons from 
thermal to relativistic energies. Coupled time-dependent diffusion equations for the electron 
distribution and the wave spectral density are solved numerically ([Miller, LaRosa, fc Moore 


19961, hereafter MLM). Nonthermal electron distributions are clearly obtained under a range 


of parameters in the accretion disk environment, which will have a direct bearing on the 
future modeling of the hard X-ray and gamma-ray emissions from GBHCs. 


THE MODEL 


2.1. Basic Model Assumptions and Parameters 


The hydrogen plasma flowing around a black hole with mass M is assumed to form 
an accretion disk (e.g., jShakura fc Sunyaev 19731) , from which most of the soft photons 
(a few keV) originate. In the inner part of the disk, we assume that the magnetic energy 
density Ub = Bq/Stt is in equipartition with the plasma thermal pressure n^ksTp, where 
k-Q is Boltzmann’s constant; Tp = Ti + Tg, and T, and (~ 50 keV) are the proton and 
electron temperatures, respectively. Ti is, unfortunately, poorly determined and is chosen 
to be 10 MeV in this study. The plasma density no is taken to be ~ l/^arR), where ctt is 
the Thomson cross section. Since Ti^ T^., it roughly corresponds to the two-temperature 
accretion disk situation. However, there have been discussions in the literature questioning 
whether a two-temperature plasma can occur at all in the accretion plasma ([Phinney 1981|; 
Rees et ah 1982|; [Begelman fc Chiueh 1988|) . We emphasize that we use Tj = 10 MeV 
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only to get a fiducial magnetic field value, and in fact, our model works even better for an 
isothermal accretion disk because electrons have a much higher thermal speed than protons 
and are preferentially accelerated (see below). 


A tenuous, quasi-spherical extended corona surrounding the hot inner disk is postulated 
(e.g., [Liang fc Price 1977| ; [Haardt &: Maraschi 199 Ij ; [Haardt, Maraschi, fc Ghisellini 1994 )- 
The coronal electron temperature is again taken to be 50 keV initially, and the coronal 
optical depth = tLc-UtR ~ 1 is varied, where ric is the coronal electron density. We further 
assume that the magnetic held in the corona is the same as that in the disk Bq, so that 
the corona is magnetically dominated and the plasma (3 = ncksTp/U-Q is < 1. Notice 
that f3 = Tc with these conditions. Adopting the black hole mass M to be IOMq and the 
size R of the system to be ~ 30GM/c^, within which most of high energy emissions are 
produced, we hnd that the coronal magnetic held Bq and dimensionless Alfven speed va/c 
are ~ 3.7 x 10® G and Bq/^AT^ njripC^ = 0.15/3“^/^, respectively. The (3 <1 condition also 
implies that the proton thermal speed is always less than va- 


Here, we make a key assumption that a fraction of the total available energy goes into 
generating MHD turbulence. This wave turbulence consists of fast mode and shear Alfven 
waves (the slow mode will be heavily Landau damped and will probably not be excited), 
but it is just the fast mode waves that are relevant for electron acceleration. Fast mode 
waves propagating obliquely with respect to the ambient magnetic held have a parallel 
magnetic held component, which can couple strongly (or resonate) with a particle when the 
parallel phase speed of the wave uj/k\\ is about equal to the parallel component of particle 
velocity n||; i.e., when vn = va/Vj where we have used the dispersion relation for the fast 
mode waves and rj is the cosine of the wave propagation angle. This process is referred to 
as transit-time damping or magnetic Landau damping (e.g., jStix 199^ ; [Lee fc Volk 1975| ; 
Achterberg 198 Ij) , and a broad-band wave spectrum leads to both rapid wave damping 
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and particle acceleration (MLM). This process is essentially the resonant form of Fermi 
acceleration ([Fermi 1949|) . Particle interactions with the parallel magnetic field variations 
(cansed by the waves) can be viewed as either head-on (gaining energy) or trailing (losing 
energy) “collisions”. Head-on collisions occnr more often than trailing ones, so that net 
acceleration resnlts. The resonant condition implies that the acceleration threshold is v\. 
For the coronal plasma, wa is comparable to the electron thermal speed but greater than 
the proton thermal speed if (3 is less than 1, so that electrons are preferentially accelerated 


2.2. The Electron and Wave Diffusion Equations and Energy Transfer Rates 


To simplify the calculations, we consider an isotropic fully-ionized hydrogen plasma 
permeated by a homogeneous background magnetic field Bq. At some large scale Ainj, 
an unspecihed mechanism generates fast mode waves. The evolution of the electron 
distribution N{E) is given by the Fokker-Planck equation 


dN 


A 



1 02 
^ 2 ^ 


[{D + D,)N] 


( 1 ) 


where E is the kinetic energy. Here, we have neglected the escape and (possible) 
pair production. Those coefficients associated with wave-particle interactions are the 
systematic mean acceleration rate {dE/dt)^^^ = p~‘^d[p^vD{p)]/dp and the diffusion 
coefficient D = 2v‘^D{p), where v and p are the electron speed and momentum, and D{p) 
is the momentum diffusion coefficient (given below). The other convection term {dE/dt)\oss 
represents the sum of electron energy change rates from inverse Compton scattering and 
synchrotron (ICS) losses {dE/dt)ics and e-e Coulomb collisions {dE/dt)^ (jPermer fc Liang 
1988 ), which also give rise to diffusion (jSpitzer 1962| ; Dermer fc Liang 1988|) . These are 


the most important processes for electrons in our parameter regime. The processes that are 
neglected include: the electron-proton Coulomb interaction since it is much slower than e-e 
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(when Ti < 100 MeV); the diffusion due to Compton scatterings since the electron recoil is 
typically small and wave-particle diffusion will dominate at high energies; the energy gain 
due to the convergence of the flow since the accretion time is much longer than dynamic 
timescale (~ R/c). The momentum diffusion coefficient D{p) for transit-time damping is 
(MLM) 

D{p) = {m^cf ^ c{k) C (-) F(/io) (2) 

lo \ c / yrrieCj \vj 


where rrie is the electron mass, C = Ut/U-q, hr is the fast mode wave energy density, {k) 
is the mean wavenumber of the fast mode waves, po = v^/v and is basically an 

efficiency factor, which equals zero when v < (i-e., when resonance is impossible) and 

equals —5/4 — (1 -|- 2pg) In po + Ro + otherwise. 


The evolution of the isotropic wave spectral density Wt can be approximated by (see 
Zhou fc Matthaeus 1990|) 


OWt 


A 

dk 




ylTr + Qw<5(fc — ko) , 


( 3 ) 


where we have included a term 7 (/c) for the wave damping by the electrons and a term 
Qw for the injection. The wave damping rate 7 (fc) can be obtained from the relation 
dk'y{k)WT{k) = dE N (E){dE/dt) (cf. equation (|)). Since {dE/dt)acc (k), this 
implies that large scale waves (small k) will suffer little loss, and that rapid wave dissipation 
only occurs when k is sufficiently large. 

At steady state, the volumetric wave energy injection rate Qw [ergs cm“^ s“^] at 
ko = 27r/Ainj must equal to the rate at which energy is transferred to smaller scales 
\/2vAkoC^^^UB (MLM). This energy flow is eventually dissipated at higher wavenumbers 
by electrons, which in turn produce X-ray to gamma-ray emissions. This implies that 
the volumetric gamma-ray energy production rate L^/{A7r/3R^), where Ly is the 







gamma-ray luminosity, must be smaller or equal (steady state) to Qw- This gives 


^ ^ Va R ^ 3/2 

T-r ~ ^ \ S ? 

Ub c Ai 


(4) 


'mj 


where = L^/{A'kR‘^c) is the energy density of gamma-ray photons. Letting Xinj/R ~ 0.1, 
we obtain U^/U-q ^ . This implies that in order to get U^/U-q ~ 0.1 as suggested 

by the 60-1000 keV luminosity from Cyg X-1 ( Phlips et al. 199^ , C, must be ~ 0.08/9^/^, 
well within the weak turbulence limit. 


Next, we look at the electron energy gain and loss rates. The mean acceleration rate of 
electrons by the waves is given as (MLM) 



where G(/io) = T"(po) + (4/io ho ~ ho + 1)/(47^) when po < 1 and equals 0 otherwise. The 
product (/c)C has to be obtained from the simulations (see next section). The ICS losses 
can be written as {d'y/dt)ics = —(4/3)(0p/rdyn)(l + U^Yl/U-q){v■, where Tdyn = R/c 
and 0p = kBTp/rrieC^. The soft photon energy density f/ph is hxed to be the same as U-q in 
this study. Note that this loss rate may be an overestimate for mildly relativistic particles, 
thus more careful treatment of the losses will help the acceleration. 

When the Coulomb loss timescale is longer than both the acceleration and the ICS 
cooling timescales, we can dehne a critical energy 7 c at which acceleration balances radiative 
cooling, (d'y/dt) = |(d7/dt)ics|, which gives 

. ( 6 ) 

where we have utilized the fact that when p/{mec) ~ 7 c 3> 1, G{pq)/tc is almost constant 
(ranging from 0.76-1.1) for 0.1 < Tc < 1. In order to get substantial acceleration—say 
7 c ~ 10—it must be that (A;) C is > 8 x 10“"^ (cm“^). This offers a direct test of our 
simulations. 
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3. Results 

We solve equations (J^) and (^) using the Crank-Nicholson method (MLM), and 
concentrate on the time evolution of particles and waves from the start of wave 
injection until the steady state is reached, during which waves are constantly injected at 
fco = 27r/ (O.li?). This period turns out to be always less than or comparable to the dynamic 
timescale Tdyn ~ 1.5 x 10“^ sec. This validates our assumption of neglecting escape. 

The Coulomb loss is hxed at the rate that is corresponding to the initial Maxwellian. 
Even though this treatment is not self-consistent, our particle acceleration results should 
not be affected much since Coulomb loss plays a negligible role compared to the ICS 
cooling at relativistic energies. The Kolmogorov phenomenology for the wave evolution is 
assumed so that ITr = with Wq cx {Qw . The injection rate Qw is 

chosen as ~ 2.8 x lO^^iyA/c) ergs cm“^ s“^, corresponding to C = 0.2 at steady state, and 
Ut = /q°° Wok~^^^dk should be a constant for all the cases considered here; there are all 
conhrmed by the simulations to within the numerical error. 

Figure 1 summarizes the time evolution of the particle density distribution N{E) 
(upper panels) as a function of kinetic energy E, and the corresponding wave spectral 
density ITt (lower panels) as a function of wavenumber k. Three different densities are 
considered. Each plot has 16 curves in it, corresponding to 5 evenly-spaced time intervals 
in each of three periods, t = 0-0.0.04-0.16rdy„, and O.lG-O.Gr^j^^, respectively. The 
particle distributions soften hrst (e.g., curves 1 — 8 in Tc = 0.5 case), due to the fact that 
acceleration is very inefficient initially since waves have not fully cascaded (i.e., small (k) as 
shown in lower panel), and Coulomb and ICS losses dominate at high energies. As waves 
cascade over the inertial range, (k) quickly grows to a level that acceleration overcomes 
both Coulomb and ICS losses. Electrons are then energized out of the thermal background 
and the nonthermal hard tail forms. This is indicated, for example, by curves 6-15 in the 
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Tc = 0.1 case. After that, both the particle and wave spectra gradually reach steady states. 

The Tc = 0.1 case clearly illustrates several other points. From simulations, {k)( ~ 10“^ 
(cm“^), thus 7 c ~ 10 using equation (|^), which is in perfect agreement with curve 15. 
Furthermore, the nonthermal tails start to develop only at E/rrieC^ ~ 0.13 (corresponding 
to va/c = 0.46), complying with the acceleration threshold. The threshold energies for 
Tc = 1 and 0.5 are buried in the thermal distributions. 

That the particle’s high energy cutoff gets larger as Tc decreases can be understood 
from equation (|^). As Tc decreases, fewer particles are available to absorb the wave energy, 
which results in larger (k). Smaller Tc also reduces the number of particles with v > va, as 
evident from the upper panels in Figure 1. 

Figure 2 shows the fraction of electrons with E > 511 keV out of the total 
electron population as a function of coronal optical depth Tc = 0.1-1, both in energy 
content /e 5 ii = [I^iEN{E)dE] / [J^ EN{E)dE] (upper panel), and in number 
/v 5 ii = [/ 5 T 1 / [ItT E^iE)dE] (lower panel). The horizontal dashed line indicates 

the initial values for a 50 keV Maxwellian. Again, we can see that a signihcant fraction of 
particles are accelerated into a nonthermal population for < 1. 


4. Conclusions and Discussions 

We have studied particle acceleration in galactic black hole accretion disk coronae 
via interactions between electrons and fast mode waves—specihcally, via the transit-time 
damping process. Including Coulomb collisions, inverse Compton scattering and synchrotron 
losses, we show that particles with speeds higher than the Alfven speed can be accelerated 
out of the thermal background, and we obtain steady state particle distributions composed 
of a Maxwellian plus a nonthermal high energy tail extending into several tens of MeV. 
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Detailed radiation modeling will be presented in a forthcoming work and we expect that 
the Maxwellian and the nonthermal tail of the particle distribution will be responsible for 
the power-law spectra in the tens of keV and the high energy gamma-rays observed from 
several GBHCs such as Cyg X-1 and GRO J0422, respectively. 

The generation of plasma wave turbulence in accretion disk environments is a fairly 
unexplored topic, but it is reasonable to suppose that the fast mode waves will be excited 
since it is an intrinsic long-wavelength mode of a magnetized plasma. We emphasize that 
the coronal plasma (3 must be < 1 for electrons to get most of the wave energy, otherwise 
proton acceleration becomes possible, reducing the energy flow to the electrons. 

The particle acceleration mechanism discussed here also has direct implications on the 
high energy radiation from accretion disks in AGNs, notably Seyfert galaxies. Preliminary 
analyses have indicated that most of our results are insensitive to the size of the system; 
thus, we expect > MeV emissions are also being produced in Seyferts as well, though high 
quality spectra above 200 keV are clearly needed to firmly settle this issue. 

We thank the referee for many constructive suggestions. Gomments from Drs. Edison 
Liang and Ghuck Dermer are appreciated. H.L. gratefully acknowledges the support of the 
Director’s Postdoc Fellowship at LANL. Part of the work was completed during J.A.M.’s 
visit at LANL, also supported by the same fellowship. 



- 12 - 

REFERENCES 

Achterberg, A. 1981, A&A, 97, 259 
Begelman, M.C. & Chiueh, T. 1988, ApJ, 332, 872 
Dermer, C. D., & Liang, E. P. 1989, ApJ, 339, 512 
Dermer, C.D., Miller, J.A., & Li, H. 1996, ApJ, 456, 106 
Fermi, E. 1949, Pliys. Rev. 75, 1169 

Gary, S.P. 1993, Theory of Space Plasma Microinstabilities (Cambridge: Cambridge 
University Press) 

Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 

Haardt, F., Maraschi, L., & Ghisellini, G. 1994, ApJ, 432, L95 

Harmon, B. A. et ah 1994, in The Second Compton Symposium, ed. C. Fichtel, N. Gehrels, 
& J. Norris (New York: AIP Conf. Proc. No. 304), 210 

Lee, M. A., & Volk, H. J., ApJ, 198, 485 

Li, H., Kusunose, M., & Liang, E. P. 1996, ApJ, 460, L29 

Liang, E. P. 1993, in The first Compton Symposium, ed. M. Friedlander, N. Gehrels & D. 
Macomb (New York: AIP Conf. Proc. No. 280), 396 

Liang, E.P. & Price, R.H. 1977, ApJ, 218, 247 

Ling, J. C. et ah 1996, ApJ, to be submitted 

McConnell, M. et ah 1996, ApJ, to be submitted 

Miller, J.A., LaRosa, T.N., & Moore, R. L. 1996, ApJ, 461, 445 (MLM) 



- 13 - 

Phinney, E.S. 1981, in Plasma Astrophysics, ed. T.D. Guyenne & G. Levy (ESA SP-161), 
337 

Phlips, B. F. et al. 1996, ApJ, 465, 907 
Rees, M.J. et al. 1982, Nature, 295, 17 

Roberts, D. A., & Wiltberger, M. J. 1995, J. Geophys. Res., 100, 3405 
Shakura, N.I. & Sunyaev, R.A. 1973, A&A, 24, 337 
Shapiro, S.L., Lightman, A.P., & Eardley, D.N. 1976, ApJ, 204, 187 
Skibo, J.G. & Dermer, G.D. 1995, ApJ, 455, L25 

Spitzer, L. 1962, Physics of Fully Ionized Gases (New York; Interscience) 

Stix, T. H. 1992, Waves in Plasmas (New York: AIP), pp. 273-292 
Sunyaev, R., & Titarchuk, L. 1980, A&A, 86, 121 
Titarchuk, L. 1994, ApJ, 434, 570 
van Dijk, R. et al. 1995, A&A, 296, L33 

Verma, M. K. 1994, Magnetohydrodynamic Turbulence Models of Solar Wind Evolution, 
Ph.D. Dissertation, Univ. of Maryland at Gollege Park 

Zhou, Y., & Matthaeus, W. H. 1990, J. Geophys. Res., 95, 14881 


This manuscript was prepared with the AAS DTeX macros v4.0. 




Fig. 1.— The time evolution of the electron density distribution as a function of kinetic 
energy E (upper panels) and wave spectral density as a function of wavenumber k (lower 
panels). Three coronal plasma densities are considered. There are 16 curves in each 
plot and the numbers indicate their time-sequence, spanning from f = 0 — 0.6Tdyn, where 
Tdyn ~ 1.5 X 10“^ sec. The initial 50 keV Maxwellians are shown as the thick-dashed curves 
(wave spectrum is zero initially) and the thick-solid curves show the steady state particle 
and wave distributions. As wave cascade progresses, the mean wavenumber becomes large 
enough that efficient acceleration results, as indicated by the development of the hard tails 
beyond the Maxwellian in particle distributions. 


Fig. 2.— The fraction of electrons with E >hll keV out of the total electron population at 
steady state for the coronal optical depth Tc ranging from 0.1 to 1, in both the energy content 
/e 5 ii (upper panel) and the number density /n 5 ii (lower panel). The dashed lines indicate 
the values for a 50 keV Maxwellian. When Tc is high, electrons are mostly nonrelativistic 
due to cooling, but a good fraction of electrons becomes relativistic and nonthermal from 
the wave acceleration when Tc < 0.5. 










